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Abstract — In wireless networks, energy consumed for commu- 
nication includes both the transmission and the processing energy. 
In this paper, point-to-point communication over a fading channel 
with an energy harvesting transmitter is studied considering 
jointly the energy costs of transmission and processing. Under the 
assumption of known energy arrival and fading profiles, optimal 
transmission policy for throughput maximization is investigated. 
Assuming that the transmitter has sufficient amount of data in its 
buffer at the beginning of the transmission period, the average 
throughput by a given deadline is maximized. Furthermore, a 
"directional glue pouring algorithm" that computes the optimal 
transmission policy is described. 

I. Introduction 

Battery size is one of the main bottlenecks on the network 
lifetime in wireless sensor networks. Replacing batteries may 
be expensive or inconvenient for nodes that are deployed 
in remote locations. In recent years, energy harvesting (EH) 
has become a viable solution to operate wireless sensor 
nodes in a self -powered fashion for extended periods of time. 
However, due to the physical and technological limitations 
of EH devices, harvested energy is typically low. Therefore, 
management of harvested energy is essential. 

In wireless systems, energy consumption for communication 
has two components: transmission energy used by the power 
amplifier and the processing energy cost (lj. Depending on the 
range of communication or the complexity of the processing 
circuitry, either of these components can be the dominating 
factor. For an energy limited system, it is known that increas- 
ing the transmission time and lowering the transmission power 
as much as possible is throughput optimal when the processing 
energy cost is ignored @. On the other hand, it is shown in 
[3| that when the processing energy cost is taken into account, 
the optimal transmission scheme becomes bursty, as increasing 
the transmission time means increasing the energy spent for 
processing. 

While most of the previous work on EH communication 
systems focus mainly on the maximization of throughput while 
ignoring the processing energy cost, in this paper we consider 
the power used by both the power amplifier and the processing 
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circuitry, and optimize the data transmission schedule over 
a fading channel. We consider a constant processing energy 
cost per unit time whenever the transmitter is operating. We 
assume that the transmitter has a finite capacity battery and 
that a sufficient amount of data is already available at the 
transmitter's data buffer before transmission starts. We focus 
on offline optimization, that is, the energy arrival instants 
and amounts as well as channel gains until the transmission 
deadline are known in advance by the transmitter. Note that the 
noncausal knowledge of energy arrivals can model accurately 
systems with predictable energy arrivals 0), or provide an 
upper bound on the performance for the case of unpredictable 
energy arrivals. Our goal is to identify an optimal transmission 
policy that maximizes the total transmitted data by a given 
deadline T subject to the energy causality constraint. 

Recently, offline transmission policies have attracted sig- 
nificant interest ll5ll- lfT4ll . Optimal transmission policies which 
account only for the power amplifier are studied in l5ll- lfT2ll . 
Yang and Ulukus Q investigate a single link EH system 
over a constant channel for given energy and data arrival 
profiles, and provide an algorithm which computes the optimal 
transmission policy. Other communication scenarios with EH 
nodes that have been studied include single link fading channel 
J6), multiple access channel J7J, interference channel and 
two-hop networks l9l- |[T0l . Battery imperfections for a single 
link system is investigated in 1111 and lfl2l . While ifTTTl 
studies battery leakage and reduction in battery capacity over 
time, 11211 investigates finite size battery. Optimal transmission 
policies over a constant channel while accounting for both 
transmission and processing energy costs are studied in lfl3l 
and HI. 

II. Preliminaries 

A. System Model 

We consider an EH communication system that harvests 
energy in packets of finite amount at time instants i e .o = < 
te.i < ■ ■ ■ < te.n < T such that the packet harvested at t = t e ^ 
has energy E^. Harvested energy is stored in a finite battery of 
capacity E max before it is used for transmission. Therefore, 
without loss of generality, we can assume that each energy 
packet can have at most E max amount of energy. We assume 
that there is no energy loss in storing and retrieving energy 
from the battery. We also assume that the real valued channel 
gain h(t) changes at time instants tf t o = < tf t i < ■ ■ ■ < 



tf.m < T, and remains constant in between. The channel is 
modeled as having additive white Gaussian noise with unit 
variance. Without loss of generality instantaneous transmission 
rate is given by Shannon capacity \ log(l + h{t)p(t)), where 
p{£) is the transmission power at time t. We assume that the 
transmitter is able to change its data rate instantaneously by 
changing the transmission power, p(t). 

We can combine all energy arrivals and changes in the 
channel gain in a single time series to = < ti < ■ • • < 
<at_i < T. This can be achieved by allowing zero energy 
arrivals at some ti's, or the channel gain to remain constant 
across some of the intervals. For consistency in notation, we 
assume an energy arrival of En = at t = T, The time 
interval between two consecutive events is called an epoch, 
and Ti = ti — denotes the duration of the i'th epoch. The 
channel gain for epoch i is denoted by hi. We are interested 
in offline optimization, that is, we assume that the transmitter 
knows all the energy arrival instants and amounts as well as 
the channel gains for the period < t < T in advance at 
t = 0. 

We assume that the transmitter consumes energy only when 
it is transmitting, and the processing energy cost is e joules 
per unit time, independent of the transmission power, p(t). 
We ignore the cost of switching the transmitter 'on' and 'off, 
and assume that no information is conveyed by the state of 
the transmitter as in 0. 

A transmission policy refers to a power allocation function 
p{t) for < t < T. A feasible transmission policy should 
satisfy the energy causality constraint: 



E(t)< 



E 

i:0<ti<t 



Ei, We[0,T], 



(1) 



where E{t) is the total consumed energy by transmission pol- 
icy pit) up to time t, i.e., E{t) = J* (p(r) + e ■ l {p(r)>0} )dT. 
In addition, battery overflows lead to a suboptimal transmis- 
sion policy because we can always increase the throughput 
by increasing the transmit power such that there is no battery 
overflow. Therefore, an optimal transmission policy must also 
satisfy the following constraint: 



]T Ei - E(t) < E max , Vie[0,T] 



(2) 



i:0<ti<t 



Assuming that the transmitter has sufficient data in its data 
buffer at time t = 0, our goal is to maximize the throughput 
under the above constraints by deadline T for given energy 
arrival and fading profile. 

B. Single Energy Arrival and Fading Level 

It is well known that for a fixed energy budget and 
no processing energy, increasing the transmission duration 
strictly increases the throughput if the rate-power function is 
non-negative, strictly concave and monotonically increasing, 
properties that are satisfied by most of the practical rate- 
power functions J2). However, if the processing energy is 
not negligible, increasing the transmission duration does not 
increase the total amount of transmitted data after a certain 



point in time since the processing energy consumption starts 
to dominate the consumed energy. It is shown in [;3] that there 
is an optimal transmission duration and power level which 
depend only on the processing energy cost e and the channel 
gain. Note that constant power level is optimal due to the 
concavity of the rate-power function. 

For a single energy packet arrival E at time t = and 
a static channel state h, we first assume that there is no 
transmission deadline. Denoting the total transmission duration 
by 9, maximum throughput is given by the solution of the 
following optimization problem: 



max — log(l + /it>), 

0,v.0(v+<i)<E 2 



(3) 



where v is the transmission power. Setting 9 = -^j- in ([3]) 
and differentiating with respect to v, the optimal transmission 
power v* should satisfy 

1 1 



-log(l + hv*). 



(4) 



Above equation has only one solution for the optimal power 
level v* which increases as the channel gain h decrease^ 
Moreover, v* does not depend on the available energy E. 
When there is a transmission deadline T, if T > v * E +e , then the 
above solution is still optimal. On the other hand, if T < v F +e , 
transmitting at power v* cannot be optimal because some 
energy would remain in the battery at time T. In this case, 
we can increase the throughput by increasing the transmission 
power so that all the available energy is consumed by time T, 
and the optimal transmission power is given by ^ — e. 

C. Related Work 

1) Glue Pouring: For a battery limited node with pro- 
cessing energy cost, when there are multiple fading levels, 
optimal transmission policy is different from the well-known 
waterfilling solution and called "glue pouring" 0. For ease 
of exposure, we describe glue pouring for two fading levels, 
single energy arrival and no deadline. Using differential power 
allocation (see for details), for single energy arrival E 
and fading states hi > h 2 with durations n, T2, respectively, 
the optimal transmission policy is summarized below. In the 
following, 9i and 92 are transmission durations for epochs 
with fading levels h\ and /12, and v* and v 2 are the solutions 
of for channel gains hi and h-2, respectively. 

• If E < Ti ivl + e), then optimal transmission policy is 
@i = an d ©2 = with power levels v\ and 0, 
respectively. 

• If Ti(u* + e) < E < Ti(v2 + — -fir + e), then optimal 
transmission policy is 9i = r\ and 92 = with power 
levels — — e and 0, respectively. 

• lfTiiv* + 1 L-± + e)<E<Tiiv* 2 + ±-± + e) + 
T 2( v 2 + e )' men optimal transmission policy is 9i = ti 
and 99 = , 2 1 



with power levels v 2 



■£7 — ^7 and u|, respectively. 



'This follows from @ by taking the derivative of v* with respect to h. 



. If n(v% + 7^ _ 7TT + e ) + T 2( u 2 + e) < -E> then optimal 
transmission policy is the usual waterfilling. 

2) Directional Waterfilling: The directional waterfilling al- 
gorithm, introduced in |6| for an EH fading communication 
system with no processing energy cost, is an adaptation of the 
classical waterfilling algorithm to the EH model where the 
energy becomes available over time. Due to energy causality 
harvested energy Ei can only be allocated to epochs j > i\ 
and, due to the battery constraint, the amount of energy that 
can be transferred to epoch j is limited by E max — Ej-i- 

III. Throughput Maximization 

A. Problem Formulation and Solution 

In this section we study the throughput maximization prob- 
lem with multiple energy arrivals and fading levels. It is 
possible to show that within each epoch, when the transmitter 
is 'on', constant power transmission is optimal J5), so we 
denote the nonnegative power level within epoch i as pi with 
duration 8;, < 0, < Tj. Then, the throughput optimization 
problem can be stated as follows: 



JV 



e, 



max — - log(l + hiPi) (5a) 

i 

s.t. 0<J2(E j -i-Q j (p j + e)),i=l,...,N, (5b) 



j'=i 



53 E i-l ~ E 6 jfe + £ ) ^ E max,i = 1, N, 



0=1 



= 1 



(5c) 

< &i < n, and < p,, i = 1, JV. (5d) 

Note that this is not a convex optimization problem because 
the constraints in d5bT>-d5cb are not convex. Therefore, we will 
reformulate this problem by defining a new variable ctj = 
QiPi, which denotes the total consumed energy by the power 
amplifier within epoch i. Then, the optimization problem in 
(0 can be written in terms of 9; and ati as follows: 



hidii 

™% 2.T log l 1 + ~0 

i—i 



= 1 

i+1 

E^-E 

= 1 = 1 



(6a) 

i = l,...,N, (6b) 



convex feasible set, therefore, the optimization problem in © 
is a convex optimization problem. 

The Lagrangian of (O with A; > 0, /i; > 0, ji > 0, v j > 
and o~i > can be written as: 



JV 



£ = E T log 1 



i=l 



(7) 



f=l \j=l 

JV / i+1 i 

" E ^ E - E (°i + e9 j) - ^" 

i=l \j=l j=l 

JV JV JV 

~ y^7»(Q» - t j) + e + e 

i—1 i—1 i—1 

Corresponding complementary slackness conditions are 

A, ( a o + eQ o - E o-i)j = 0, Vz (8) 

£ ^j-i - E ^ + e0 ^ ~ Ema A = °> Vi (9) 

7i(©i-Ti)=0, ^6j=0 and ^a* =0, V*. (10) 
Taking derivatives with respect to a, and 6^, we obtain 

JV 



dC 



Qihi 



don 2(0, + /ijCKi) 
<9£ 1 



JV 



(12) 



We consider KKT conditions together with the complementary 
slackness conditions in d8l-([T0li. 

• If 0* = 0, then a* = and no power is allocated to 

epoch i, i.e., p* = 0. 
. If < 0* < n and < a*, i.e., 7i = 0, u t = 0, cr 4 = 0, 
then from ( fTTT i and ( fT2l , we obtain 



s.t. < 53 (-Ey-i - a o ■■ ~ e0 i)i 

(Oj- + e0j) < £7 TOOX , « = 1, —, AT, 



10g 1 0* ) 0? + fca? ' 



(13) 



(6c) 

< 0j < ri, and < q„ i = 1, JV. (6d) 

With this reformulation, concavity of < f6ab can be argued 
from the fact that the function ^Mog(l + is the per- 

spective of the strictly concave function 4 log(l+/ijOi). Since 
perspective operation preserves concavity, ^ log(l + is 
also concave ifTBI . The linear constraints in (l6bT>-(f6db define a 



which is equivalent to (01i when we replace a* with 0*p;. 
Therefore, as argued in Section III-B1 ( fl3l l has a unique 
solution which depends only on hi and e. We denote this 
unique power level as p* = v* and note that p* does 
not depend on the Lagrange multipliers Xj and fj,j, for 
j = i, N, either. 

If 0* = ^ and < a*, i.e., j t > 0, v. t = 0, a l = 0, then 
from ( fTTI ). we get 



2(0? + M?) 



(14) 



0=1 



Similarly, from ( TTZb . we get 



hid* 



hiQ* 



9* y 2(9* + M?) 

JV 

+e^(A i -^)+7 i . (15) 

Since 7^ > 0, from ( fT4t and dT5T >. we can obtain the 
following inequality: 



log 1 



> 



(16) 



e* ) 0* + h ia * 

Comparing with ( TT3b we conclude that ( [ToT l is satisfied 
only if p* > u*. Then, by using (fl4l i and replacing a* 
with 9*p*, we can compute the optimal power level p* 

as 



1 



2 EL (Ai-Af,-) 



(17) 



Notice that Xj and /ij cannot be positive simultaneously. 
Using the complementary slackness conditions in ([H) and 
(0, we have A, > and /i, = whenever the battery 
of the transmitter depletes. Therefore, from (fl7l l. we can 
argue that + ^— j- > p* + This means that 
whenever the sum of the inverse channel gain and the 
optimal power level increases from one epoch to the 
next, the battery must be empty. When /Xj > and 
Aj = 0, the battery is full. Therefore, from ( fTTI i. we 
can argue that p* +1 + ^— j- < p* + As a result, the 
sum of the inverse channel gain and the optimal power 
level decreases from one epoch to the next, whenever the 
battery is full. Moreover, since depleting all the harvested 
energy by the deadline is optimal IT2l . Ajv > and 

MTV = 0. 

Remark 3.1: Optimization problem in © may have multi- 
ple solutions. Consider a channel with multiple epochs having 
the same channel gain hi and a corresponding optimal trans- 
mission policy with < 9* < Tj. As argued from KKT 
conditions, p* = v* must be satisfied for these epochs. Then 
the corresponding optimal values for = hip* must also 
be the same. Without loss of optimality, we can find another 
optimal transmission policy by transferring some of the energy 
between those epochs such that the optimal hip* is preserved. 
Note that if the channel gains are different (Section Ill-Bl i or 
p* > v*, i.e., 9* = n or A, > or ^ > 0, for all i, then 
there is a unique solution. 

B. Directional Backward Glue Pouring Algorithm 

We can allocate the harvested energy to epochs starting from 
the last non-zero energy packet to the first such that constraints 
in ([Hi-© (akin to directional waterfilling in Section Hl-Cl l are 
satisfied. In addition, the optimal transmission policy utilizes 
epoch i either partially, i.e. 0* < r^, with power level p* — v*, 
or fully, i.e., 0* = r,, with power level p* > v* as argued 
in Section IIII-AI Therefore, the optimal transmission policy is 
a directional glue pouring algorithm in which each harvested 



energy packet Ei is allocated to subsequent epochs using the 
glue pouring algorithm. 
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Fig. 1. Directional backward glue pouring algorithm 



Consider the example in Fig. |l(a)| Arrival times of the 
harvested energy packets are shown with thick downward 
arrows. The thin downward arrows correspond to the time 
instants when the channel gain changes with virtual energy 
arrivals (E\ = E3 = E4 = 0). Inverse of the channel gains 
are indicated by solid blocks. Optimal power levels v* for 
< 0^ < Ti are indicated with dashed horizontal lines above 
the inverse channel gain blocks such that v* corresponds to 
the distance between the dashed lines and solid blocks. Note 
that the energy consumed for processing is not shown in the 
figure; however, it can be computed from the total transmission 
duration. As argued before, the algorithm first computes the 
optimal power level for the last non-zero energy arrival i?2- 
As shown in Fig. |l(b)| the algorithm considers the harvested 
energy E2 for epochs three, four and five. It allocates E2 to the 
third and fourth epochs using glue pouring algorithm as argued 
in Section III-CI Note that only a portion of the third epoch 
is utilized with power level v$ due to glue pouring. Then, the 
algorithm considers the first non-zero energy arrival Eq and 
allocates this energy according to glue pouring algorithm as 
shown in Fig. |l(c)| Note that some of the energy is transferred 
to third and fourth epochs as argued in Section ITl-CI However, 
transferred energy is limited due to the finite battery size which 
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Fig. 3. Average throughput versus processing energy cost. 
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Fig. 2. (a) Optimal power levels for e = are P = [2.17,0.29,2.01, 
1.18,1.65] pW with durations = [0.5,3.5,1.1,1.9,3.0] sec. Total 
transmitted data is B = 2.11 nats. (b) Optimal power levels for e = 
1 pW are P = [1.99,3.48,3.05,0,1.99] pW with durations & = 
[0.36, 0.22, 1.10, 0, 1.66] sec. Total transmitted data is B = 1.39 nats. 

explains the water level difference (Section Ullt. 

IV. Numerical Results 

In this section, we provide numerical results to show the 
effect of processing energy cost on the optimal throughput. We 
consider an energy arrival profile E = [1.1,3.2,2.8,1.4,3.1] 
microjoules (fxJ), channel gains h = [0.7,0.2,0.4,0.3,0.7] x 
10 6 with epoch durations t = [0.5,3.5,1.1,1.9,3.0] sec 
until deadline T = 10 sec. We consider that energy packets 
arrive at time instants when channel gain changes. We set 
Emax = 5 piJ. The optimal offline transmission policy with 
no processing energy cost, i.e., e = is shown in Fig. 
|2(a)| As we can see from Fig. |2(a)| the transmitter utilizes 
each epoch fully. The difference in power levels is due to 
the energy causality and finite battery capacity. Using the 
same energy arrival and channel profile for e = 1 fiW, we 
obtain the transmission policy in Fig. |2(b)| As shown in the 
figure, the optimal transmission policy is bursty, and the total 
transmission energy is reduced due to the processing energy 
cost. Notice that the optimal policy allocates energy to epoch 
two even though it has the worst channel gain. This is due to 
finite capacity battery. 

The variation of the throughput with respect to e for the 
same energy and channel profile given above is shown in 
Fig. [3] As expected the optimal throughput decreases as the 
processing energy cost increases. 



V. Conclusions 

In this paper, we have studied an EH communication system 
with processing energy cost over a fading channel. Under the 
noncasual knowledge of energy packet arrivals, we have iden- 
tified the optimal transmission policy such that the transmitted 
data is maximized by a given deadline. Our solution involves 
a convex optimization formulation of the problem as well 
as an optimal 'directional glue pouring' algorithm. Finally, 
numerical results have been provided to illustrate the effect of 
processing energy cost on the throughput. 

References 

[I] S. Cui, AJ. Goldsmith, and A. Bahai, "Energy-constrained modulation 
optimization," IEEE Trans. Wireless Commun., vol. 4, no. 5, pp. 2349- 
2360, Sep. 2005. 

[2] S. Verdu, "On channel capacity per unit cost," IEEE Trans. Inf. Theory, 
vol. 36, no. 5, pp. 1019-1030, Sep. 1990. 

[3] P. Youssef-Massaad, L. Zheng, and M. Medard, "Bursty transmission and 
glue pouring: On wireless channels with overhead costs," IEEE Trans. 
Wireless Commun., vol. 7, no. 12, pp. 5188-5194, Dec. 2008. 

[4] M. Gorlatova, A. Wallwater, and G. Zussman, "Networking low-power 
energy harvesting devices: Measurements and algorithms," Proc. IEEE 
INFOCOM. Shanghai, China, Apr. 2011. 

[5] J. Yang and S. Ulukus, "Optimal packet scheduling in an energy harvest- 
ing communication system," IEEE Trans, on Communications, vol. 60, 
no. 1, pp. 220-230, Jan. 2012. 

[6] O. Ozel, K. Tutuncuoglu, J. Yang, S. Ulukus, and A. Yener, "Transmission 
with energy harvesting nodes in fading wireless channels: Optimal 
policies," IEEEJSAC, vol. 29, no. 8, pp. 1732-1743, Sep. 2011. 

[7] J. Yang and S. Ulukus, "Optimal packet scheduling in a multiple access 
channel with rechargeable nodes," Proc. ICC, Kyoto, Japan, Jun. 2011. 

[8] K. Tutuncuoglu and A. Yener, "Optimal power control for energy harvest- 
ing transmitters in an interference channel," Proc. 45th Asilomar Conf. 
Signals, Systems and Computers, Pacific Grove, CA, Nov. 2011. 

[9] D. Gunduz and B. Devillers, "Two-hop communication with energy 
harvesting," Proc. CAMSAP, San Juan, PR, Dec. 2011. 

[10] O. Orhan and E. Erkip, "Optimal transmission policies for energy 
harvesting two-hop networks," Arxiv preprint arXiv:1202.6654 Mar. 
2012. 

[II] B. Devillers and D. Gunduz, "A general framework for the optimization 
of energy harvesting communication systems with battery imperfections," 
Journal of Commun. and Netw., Spec. Issue on Energy Harvesting in 
Wireless Netw., vol. 14, no. 2, pp. 130-139, Apr. 2012. 

[12] K. Tutuncuoglu and A. Yener, "Optimum transmission policies for 
battery limited energy harvesting nodes," IEEE Trans. Wireless Commun., 
vol. 11, no. 3, pp. 1180-1189, Mar. 2012. 

[13] J. Xu and R. Zhang, "Throughput optimal policies for energy harvest- 
ing wireless transmitters with non-ideal circuit power", Arxiv preprint 
larXiv: 1204.38181 Apr. 2012. 

[14] Q. Bai, J. Li, and J. A. Nossek, "Throughput maximizing transmission 
strategy of energy harvesting nodes," Proc. IWCLD, Rennes, France, Nov. 
2011. 

[15] S. Boyd and L. Vandenberghe, Convex Optimization, Cambridge Uni- 
versity Press, 2004. 



